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WITH LOCALLY NEGATIVE DRIFT 



ANDREJ DEPPERSCHMIDT AND SOPHIA GOTZ 

Abstract. We consider two reflecting difl^usion processes {Xt)t>o with a moving reflection 
boundary given by a non-decreasing pure jump Markov process {Rt)t>0- Between the jumps 
of the reflection boundary the difi'usion part behaves as a reflecting Brownian motion with 
negative drift or as a reflecting Ornstein-Uhlenbeck process. In both cases at rate 'y{Xt — Rt) 
for some 7 > the reflection boundary jumps to a new value chosen uniformly in [Rt— , Xt]. 
Since after each jump of the reflection boundary the diffusions are reflected at a higher level 
we call the processes Brownian ratchet and Ornstein- Uhlenbeck ratchet. Such diffusion ratch- 
ets are biologically motivated by passive protein transport across membranes. The processes 
considered here are generalisations of the Brownian ratchet (without drift) studied in (Dcp- 
jjerschmidt and Pfaffclhubcr, 2010). For both processes we prove a law of large numbers, in 
particular each of the ratchets moves to inflnity at a positive speed which can be computed 
explicitly, and a central limit theorem. 



1. Introduction 

Reflecting diffusion processes constitute an important class of stochastic processes that appear 
in various applications. We consider two particular examples of a diffusion ratchet variants of 
which were introduced in (Simon ct al., 1992) and (Pcskin ct al., 1993) motivated by protein 
transport across cell membranes. Generally speaking a diffusion ratchet is a diffusion process 
reflected at a non-decreasing jump process. The name ratchet is justified by the fact that each 
jump prevents the diffusion from attaining lower values. In a sense a jump of the reflection 
boundary process can be thought of as a click of a ratchet. In (Budhiraja and Pricks, 2006) 
a diffusion ratchet (modelling a molecular motor) in which a particle moves according to a 
Brownian motion between equally spaced (deterministic) barriers is studied. The particle can 
cross such barriers from left to right but is refiected if it hits a barrier to its left. The two 
models that we consider in the present note are both generalisations of the diffusion ratchet 
studied in (Dcppcrschniidt and Pfaffclhubcr, 2010). In the model studied there a particle moves 
according to a reflecting Brownian motion and the reflection boundary jumps a rate proportional 
to the distance between the particle and the current reflection boundary. At jump times the new 
reflection boundary is chosen uniformly between the old one and the position of the particle. 

Let us now introduce the models that we consider here, then briefly explain the biological 
motivation and finally state our results. 

1.1. The models. Let {X . TV) :— {Xt, Rt)t>o be a time-homogeneous Markov process starting 
in {Xo,Ro) = {xo,0), for some xq > 0. Here X = {Xt)t>o is a diffusion process refiected at 
a non-decreasing jump process TZ = (i?t)t>o- Given that {X,Ti) is in {Xt,Rt) at time t, the 
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reflection boundary process jumps at rate ^{Xt — Rt) for some 7 > 0. If t is a jump time of the 
reflection boundary then the new position is uniformly distributed on the interval [Rt-, Xt]. By 
this dynamics, Rt < Xt for all t > 0, almost surely. In principle the above description works with 
any reflecting diffusion between the jumps. As mentioned earlier we study here two particular 
cases and to distinguish them we write (A", 7?.) and (A", 72.) for the corresponding processes. 

(I) If for /i > the process X — {Xt)t>o is a Brownian motion with negative infinitesimal 
drift — /i, unit variance (see Section 2.1) and reflection boundary process TZ = {Rt)t>o 
then we refer to the process (X,TZ) as the {■^, fi)- Brownian ratchet. 
(II) If for /i > the process X — {Xt)t>o is an Ornstein-Uhlenbeck process with infinitesimal 
drift —fix, unit variance (see Section 3.1) and reflection process TZ — {Rt)t>o then we 
refer to {X,TZ) as the {■^, fx)- Ornstein-Uhlenbeck ratchet. 
Whenever we want to stress the dependence on the parameters, we write {X'^'<'^'\^Z'^^''''^) = 
{x[^'^\R^;^'^^)t>o for the (7, At)-Brownian ratchet and {X^'^■^') ,n^^'^')) = {x[''-''\R}i''''^)t>Q for 
the (7, /i)-Ornstein-Uhlenbeck ratchet. 




Figure 1.1. A diagram of the ratcheting mechanism for protein transport on 
the left and of a reflecting diffusion with negative local drift on the right. 

1.2. Biological motivation. Inside a typical cell different proteins are involved in many pro- 
cesses. They usually need to be transported after or during the translation (production) to 
various locations at which they are required. Depending on the protein and its functions there 
are different transport mechanisms. In the present paper we focus on the passive protein trans- 
port across membranes of e.g. endoplasmic reticulum (ER) or mitochondria for which ratcheting 
models were introduced by Simon ct al. (1992) and Pcskin et al. (1993). The main idea in these 
models is that due to thermal fluctuations the protein moves, say inside and outside the ER for 
definiteness, through a nanopore in the membrane according to a diffusion; see the left part of 
Figure 1.1. Inside the ER, ratcheting molecules can bind to the protein at a certain rate. These 
ratcheting molecules are too big (in our model they are actually infinitesimally small but one 
can imagine that binding of the molecules leads to a deformation of the protein at the ratcheting 
sites) to pass through the nanopore and prevent the protein from diffusing outside the ER, i.e. the 
protein performs a reflected diffusion with jumping reflection boundary which is due to binding 
of new ratcheting molecules. In the last two decades such models have been studied extensively 
in biology, physics as well as in mathematics. For a detailed overview of the recent literature and 
for more biological motivation we refer to (Depperschmidt et al., 2012) and references therein. 

With this motivation in mind Xt (and Xt) can be interpreted as the length of the protein 
inside ER at time t and Rt (and Rt) as the distance between the "head" of the protein and 
the ratcheting molecule closest to the nanopore; see the right part of Figure 1.1. Since typically 
proteins have to be unfolded during translocation into ER, the movement inside takes place 
against a force pointing outside which explains the locally negative drift of the ratchets. 
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1.3. Results. For both, the (7, ^)-Browniaii ratchet and the (7, /i)-Ornstein-Uhlenbeck ratchet 
we prove a law of large numbers as well as a central limit theorem. Furthermore we compute 
the speed of the ratchets in terms of the Airy Ai-function in the case of Brownian ratchet and 
in terms of the Tricomi confluent hypergeometric function in the case of Ornstein-Uhlenbeck 
ratchet. 

Theorem 1.1 (LLN and CLT for the Brownian ratchet). Let {X,TZ) — {Xt, Rt)t>o be the 
{■y, Ij)- Brownian ratchet starting in (a;o,0) with a^o > 0. //7,/i > then 

T ^"(^'^) •"~2^ ^z((27)-V3,,2) '2^^ almost surely, (1.1) 

where Ai{-) is the Airy function. Furthermore in the case 7 > there is a = a(fi,y) > such 
that 

Here " =^ " denotes convergence in distribution and X is a standard Gaussian random variable. 

Note that though the result is formulated for /i > for the proof we only need to consider the 
case > 0. In the case /i = the (7, /i)-Brownian ratchet as well as the(7, /i)-Ornstein-Uhlenbeck 
ratchet reduce to the process studied in (Deppcrschmidt and Pfaffelhubcr, 2010). 

Theorem 1.2 (LLN and CLT for the Ornstein-Uhlenbeck ratchet). Assume fi > and 7 > 0. 
Let {X,TZ) =^ {Xt, Rt)t>o be the (-j, fi)- Ornstein-Uhlenbeck ratchet .starting in (a;o,0) for xq > 0. 
For X > set 



.(x)-e-W™VQ-^,l,(^ + V^.)^), (1.2) 
where U is the Tricomi confluent hypergeometric function (see (3.14) for a definition). Then 

T ^"(^'^)-=-2V;(0) KMO) ^^-^ost surely. (1.3) 

Furthermore in the case 7 > there is a — a{p,"f) > such that 

Xt - tv{fl, 7) t^oo 



for a .standard Gaussian random variable X. 



X 



Remark 1.3 (Comparison of the ratchets). In Figure 1.2 we plot the speed of both ratchets in 
the interval fi G [0, 8] in the case 7 = 1/2. The plots are based on numerical computations using 
Mathematica. In the neighbourhood of zero, here approximately in the interval (0,0.6), the 
Brownian ratchet is faster whereas outside that interval the Ornstein-Uhlenbeck ratchet has a 
higher speed. Heuristically this can be explained: If Xt — Rt ~ Xt — Rt are large and /i small 
then the drift of Xt towards Rt is smaller than that of Xt towards Rt . Then in the Brownian case 
the reflection boundary jumps on average "earlier" and "higher" than in the Ornstein-Uhlenbeck 
case. Since both Xt — Rt and Xt — Rt are "shortened" at rate proportional to their values the 
described effect is not very pronounced and the speed of both ratchets is comparable in this 
region. 

On the other hand if fx is large and Xt — Rt ^ Xt ~ Rt are close to zero then Rt has a higher 
chance to jump "earlier" and "higher" than Rt because the drift Xt towards Rt is constant and 
that of Xt towards Rt is proportional to their distance which is small in this case. 
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0.6 4 8 

Figure 1.2. Comparison of the speed of the ratchets for 7 = ^ and fi e [0, 8]. 

Note that the above heuristic arguments are similar in spirit to the following considerations in 
the case without jumping reflection boundaries (7 = 0, in that case the speed of both ratchets is 
zero). The invariant density of the reflected Brownian motion with negative drift —fj. is f{x) = 
a; > (see e.g. Harrison, 1985, p. 94) and that of the reflected Ornstein-Uhlenbeck 
process with drift — /ix and unite variance is g{x) = 1^ [ijne^^^ , a; > (this can be easily 
obtained from the invariant density of the Ornstein-Uhlenbeck process). For the expectations 
we have Jg xj(x)dx — l/{2^) and xg{x)dx — l/^/WJi. In particular, under the invariant 
distributions the expectation of the reflected Brownian motion is larger than that of the reflected 
Ornstein-Uhlenbeck process for /i < 7r/4 whereas for /i > 7r/4 the opposite inequality holds. □ 

Outline. The rest of the paper is split in two Sections in which Theorem 1.1 and Theorem 1.2 
are proved. In Section 2 we deal with the Brownian ratchet. First, in Section 2.1 we recall an 
explicit construction of the reflecting Brownian motion with drift. It will be used in Section 2.2 
to give a graphical construction of the Brownian ratchet. There we also prove a scaling property 
for the Brownian ratchet and show that the graphical construction can also be used to construct 
a coupling of Brownian ratchets with different initial conditions. Between the jump times the 
Brownian ratchet can be seen as a killed reflecting Brownian motion with drift. For that reason 
in Section 2.3 we compute the corresponding Green function and obtain several estimates on the 
moments of the killing time and the position at killing time. In Section 2.4 we study the Markov 
chain of the increments of the Brownian ratchet at jump times. We show that this Markov chain 
possesses a unique invariant distribution and compute the expectations under this distribution. 
These will be used later to compute the speed of the ratchet explicitly. In Section 2.5 we 
define a regeneration structure for the Brownian ratchet and show that the increment at these 
regeneration times have finite second moments. From that we obtain in Section 2.6 the assertion 
of Theorem 1.1. In Section -3, which has a similar structure to Section 2, we carry out the 
corresponding program for the Ornstein-Uhlenbeck ratchet. 

2. Brownian ratchet with negative local drift 

In this section we give a graphical construction of the Brownian ratchet with negative local 
drift from which we deduce a scaling property and show that the construction allows to couple 
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two Brownian ratchets so that from some almost surely finite time on they have the same spatial 
as well as temporal increments. Then we study the Markov chain of the increments of the ratchet 
at the jump times of the boundary and show that it has a unique invariant distribution, which will 
allow to compute the speed of the ratchet explicitly. For the LLN and CLT we define regeneration 
times of the ratchet and show that the increments between these times have bounded second 
moments. 

Before we start with the above schedule let us recall the definition and an explicit construction 
of the reflecting Brownian motion with drift. 

2.1. Reflecting Brownian motion with drift. Though the definition given here is valid for 
any /it S K we will assume /i > because the case < is less interesting. For more information 
on reflecting Brownian motion with drift we refer to e.g. (Harrison, 1985; Graversen and Shiryacv, 
2000; Peskir, 2006). A reflecting Brownian motion with infinitesimal drift — /x started in a; > 0, 
which we denote by RB {— ji) , is a strong Markov process with continuous paths (i.e. a diffusion 
process) associated with the infinitesimal operator acting on 

V{An := {/ e C2(M+) : /'(0+) = 0} 

as follows: 

^^/(y):-^/"(2/)-M/'(y). (2.1) 

We shall omit the superscript x and write RBM{—^) whenever the initial value is not important. 

Let us also recall from (Peskir, 2006) an explicit construction of i?BM^(— /i) that will be 
useful for our purposes. Let B = {Bt)t>o be a standard Brownian motion starting in 0. We 
define the Brownian motion with drift /i, denoted by B^^, and its running maximum, denoted by 
M^, by 

B^ = Bt+fit and M/' = max B^, for t > 0. (2.2) 

0<s<i 

Furthermore we define Z^'^ — {Z^'^)t>o by 

Zlt'"" ^ {xy Mt) - Bi;. (2.3) 
Then in (Peskir, 2006, Thm. 2.1) it is shown that 

RB]VF{-^i) = Z^'^. (2.4) 

2.2. Graphical construction of the Brownian ratchet with negative local drift. Assume 
/i > and let i?^ be as in (2.2). Furthermore let N'^ be an independent Poisson process on 
M X [0, oo) with intensity ^X^{dx^ dt) where is Lebesgue measure on . 

We define a sequence of jump times (t„)„=o,i,... and a sequence (5'^"^)„=o,i,... with 5^"^ = 
{S[^^)t>T^ as follows: 

TO := 0, (2.5) 

:= xa, Sf^^ := max{5^°\ sup {B,^}}. (2.6) 

0<s<t 

Given r^-i and 5*'"^^^ for some n>l we set 

r„ := inf{t > r„_i : n [Bf , s[''-^^] x {t} ^ 0}. (2.7) 

Furthermore we let S'l"'' be the space component of the almost surely unique element of iV n 
[B^^,5|r^^] X {t„}. For t > t„ define 

:=max{5(:'\ sup {B^}}. 

T„<S<t 
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Finally we define S {St)t>o and {X,TZ) :— {Xt,Rt)t>o by setting 

n 

St = S^"^ and Rt^Y. (^-r'' ~ for t € [r„, r„+i) 



(2.8) 



and Xt^Rt + St- for t > 0. 



Note that S is the "running maximum" process that jumps down to Poisson points that are 
between the process itself and the Brownian motion with drift. 



- 




Figure 2.3. Graphical construction of the Brownian ratchet with locally neg- 
ative drift. 



In the following lemma we verify that (X,TZ) fits the description of (7, /i)-Brownian ratchet 
given in Subsection 1.1. 

Lemma 2.1. The process {X,TZ) is ('y, 11) -Brownian ratchet started in (a:, 0). 

Proof. By construction sil > Srf, so that TZ is non-decreasing. Furthermore, St > Bt implies 
Xt > Rt for all t > 0. Between r„ and t„+i the process X is RBM^—fi) reflected at Rt starting 
at time r„ in 



R^^^ + Sr,^ - Bif^^ = J2is^;-^^ - 5(f) + - B^,^ 



i=l 
n-1 



- ) + 5(^-1) - B^^ = Rr„_^ + Sr„- - B!^^^ 



Thus, Xr„ = Xr^^ for all n > 0, and therefore the paths of X are continuous. 

Given the process up to time t„ the jump rate of TZ i.e. the rate of at which t„+i occurs is 



7(S'f — B^) = 7(Xt — Rt). Then the reflection boundary jumps to i?r„+i — Rt„ + sl^\^ — ■S'rl'+i 
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By homogeneity of the Poisson process , St^^P is uniform on [Bf^^^^, St,1^_^^]. Thus, i?r„+i is 
uniform on [i?.r„, ^t„+i] because for some U ~ U{[0, 1]) we have 

= Rr^ + (1 - U){S^l\^ - B^^J = Rr„ + (1 - C/)(X..+ , - RrJ. 

□ 

If we transform time and space in the graphical construction then we of course rescale the 
Brownian motion with drift and transform the Poisson process. There is only one such trans- 
formation that maps the Poisson process to N'^ and the Brownian motion with drift to a 
Brownian motion with another drift. 

Lemma 2.2 (Scaling property). For 7 > and fi > we have 

[^t jt>o - 7 \^j2/3t ^R-y2/3t )t>a- [^■y) 

Proof. Assume that we construct (X^^'^ ^ ^, R]''^ ' ^)t>o starting in (7^/^x0, 0) using the Pois- 
son process and the Brownian motion with drift B'^ ^'^^ . We define g : — > by 

On the one hand, rescaling space and time using g we obtain the process on the right hand side 
of (2.9). On the other hand the rescaled graphical construction leads to the process on the left 
hand side of (2.9). For that we need to verify 

g(N^)^N'^^ (2.10) 

ff((i,i?r"'')t>o) = (t,i?r),>o- (2.11) 

Equation (2.10) is clear. For (2.11) we have using the scaling property of the Brownian motion 

),>o = (7-^/^^,/^7-^/^t + 7-^/^i??)),>o 



g[it^Br-'"\>.) = (7-^/^^,7-^/^5r""\>o - {l-'/%,^-'^h + ^-^/^B?))^^^ 

(2.12) 



□ 



We now turn to the construction of a coupling of two (7, /x)-Brownian ratchets starting in 
((XjO), (a;, 0)). Let {B^)t>Q and N'^ be as before, and let x,x >Q with x >x without loss of 
generality. To construct the coupled Brownian ratchet 

((Xt,i?0,(^tJt))t>o' (2-13) 

set To = To = 0, S'o"'' = X, Sq^'^ = X and define as before the sequences (t„)„>o, (t„)„>o, (S'*^"^)„>o 
and (-S^"'')n>o- Furthermore define the corresponding processes S, S, {Xt,Rt)t>n and {Xt,Rt)t>o 
as in (2.8). Define the coupling time by 

Tcoupi = inf{t > : 5t = (2.14) 

Note that, since we use the same Brownian motion and the same Poisson process for both ratchets 
we have St = St for all t > Tcoupi- Thus, on the event {Tcoupi < 00}, there are n, n such that for 
fc > we have almost surely 

Xr^+k+l ~ -^T„+k = ^T„+fc + i ~ -'^Tn+fc) (2-15) 
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The following lemma shows that Tcoupi is almost sm'ely finite, i.e., the coupling is successful. 

Lemma 2.3 (Exponential moments of the coupling time). 
For /i > and < a < we have 

Proof. At time Tcoupi either both S and S use the same point of the Poisson process N'' or the 
Brownian motion B^^ touches the maximum of S and S which is S by assumption x < x. Thus, 
we have Tcoupi < T for 

T := inf{t >0:St^ AfJ = inf{i > : St = B^}. 

By construction St can increase only after this time T and decrease by jumping down when it 
uses points of the Poisson process. Ignoring this decrease by jumping down we obtain 

T <H^ := inf{t > : B^ = x}. 

It is well known that (see e.g. Borodin and Salmincn, 2002, p. 295) for a < At^/2 we have 

and the result follows. □ 

2.3. Green function of the killed reflected Brovifnian motion with drift. Since between 
the jumps the ratchet constructed in the previous section behaves as a killed reflected Brownian 
motion we will need in the sequel some functionals of that process, such as expected killing 
time or expected position at killing. To this end we need to compute the corresponding Green 
function. 

Let us first give a short description on how the Green function of a killed diffusion can be 
computed; for details we refer to Chapter II in (Borodin and Salmincn, 2002) or Chapter 4 in 
(Ito and McKean, 1974). 

Remark 2.4 (Green function of a reflected diffusion with killing). 

Let Y := {Y{t))t>Q be a reflected diffusion process with killing on state-space [0,oo) associated 
with infinitesimal operator A acting on 

V{A) {/eC,2(M+):/'(0+)=0} 

as follows 

Afiy) lf"{y) + b{x)r{y) - c{x)f{x). (2.16) 

We consider in this paper the cases b{x) = — or h{x) = —fix and c{x) = jx. Since in both cases 
the killing time is almost surely finite the resulting diffusions are transient. The speed and the 
killing measures of Y are given by (see e.g. Borodin and Salmincn, 2002, p. 17) 

m{dx) ^ m{x) dx 2e^^'''> dx and k{dx) = k{x) dx := 2c{x)e^'^''^ dx, (2.17) 

where B{x) := J^^ 2b{y) dy. Let p(-; •, •) denote the transition density of Y with respect to the 
speed measure. Then the Green function of Y is defined by 

G{x,y) / p{t;x,y)dt. 
Jo 

In the transient regular case (the latter means here that every point in [0, oo) can be reached 
with positive probability starting from any other point) the Green function of Y is positive and 
finite. It is obtained in terms of two independent solutions (p and -0 of the differential equation 
Af ~ that are both unique up to a constant factor and satisfy the following conditions 



DIFFUSION RATCHETS WITH LOCALLY NEGATIVE DRIFT 



9 



(i) (j) is positive and strictly decreasing with 4i{x) — ^ as x — > oo, 

(ii) ijj is positive and strictly increasing, 

(iii) i/j'{0+) = (this condition is for the reflecting boundary). 

The Wronskian, defined by w^ip^cf)) := ^'{x)(j){x) — ip{x)'ip'{x) is independent of x. Thus, the 
functions (f) and can be chosen so that their Wronskian equals one. Then, the Green function 
of Y is given by 

G(x,y) = |^^")^^(y) : 0<y<x, ^^^^^ 
\^tp(x)(l){y) : 0<x<y. 

In our computations we will in principle not need the exact expressions for (jj, ip and G. In 
particular in the case of the Ornstein-Uhlenbeck process, where these solutions depend in a 
complicated manner on the model parameters, we will not compute the function -i/j explicitly. 
Asymptotic bounds at infinity will suffice for our purposes. 

In the following remark we collect some properties of the Airy functions that will be needed in 
the sequel. For further properties we refer to (Abramowitz and Stegun, 1992) (cf. also Remark 5.2 
in (Dcppcrschmidt and Pfaffclhubcr, 2010)). 

Remark 2.5 (Airy functions). The Airy functions Ai and Bi are two linearly independent solu- 
tions of the differential equation 

u"{x) - xu{x) ^0. (2.19) 

We will only need the properties of the Airy functions on [0,oo). On that domain the functions 
are positive, Ai is decreasing with Ai{x) q and Bi is increasing with Biix) ^^°°) qo. 

The Wronskian is independent of x and is given by 

w{Ai, Bi) = Bi'{{))Ai{Q) - Ai'{Q)Bi{Q) = -. (2.20) 

TT 



The integral of Ai on IR+ is 



1 

Ai{u)du=-. (2.21) 
o 



We will also need the function 

Gi{x) -.^ Ai{x) Bi{y) dy + Bi{x) Ai{y)dy. (2.22) 
Ja Jx 

For fixed /i > and C G M we define a function M by 

M{x) 7r|Ai(^2 + x) ^ {Bi{i? + y) + CAiip"^ + y)) dy 

+ {Bi{^i^ +x)+ CAiiiT' + x)) I Ai{fi^ + y) dy\. 
Using (2.21), (2.22) and positivity of Ai and Bi on [0, oo) we have for x > 



(2.23) 



M{x) < TT ( Gi{p^ +x) + ^^i(/i^ + (2.24) 



Since the functions Gi and Ai are bounded on [0, oo) we may define 



G* ■.^G*{C,ii) -.^TT-mayilGiii? + x)+^^Ai{ii^ +x)\. (2.25) 

□ 
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Definition 2.6 (Killed reflecting Brownian motion with drift). Let Z := {Zt)t>o denote a 
reflecting Brownian motion with drift —fi starting in a; > (see (2.3)) under the law P^; and let 
Ea; denote the corresponding expectation. Furthermore using an exponentially distributed rate 
1 random variable ^ independent of Z we define the killing time by 

J. 



m{{t > : 7 / Zsds>^}. 
Jo 



Then reflecting Brownian motion with infinitesimal drift — /i killed at rate "fZ is defined as the 
process Z^ := {Zf)t>o with = Zt for t £ [0, t) and Zf = A for t > r for some A ^ M, often 
referred to as the cemetery state. The infinitesimal operator A^''' corresponding to Z^ acts on 
functions / : [0, oo) — )■ M satisfying /'(0+) = as follows 

A^^-^ix) = \r{x) ~ t^f'ix) ~ ^xfix). (2.26) 

In view of the scaling property it is enough to prove the results for the (7, /L()-Brownian ratchet 
in a particular case. In what follows we assume 7=5 and fix /i > 0. In this case the speed and 
killing measure corresponding to Z^ are given by 

m{dx) = m{x) dx := 2e^^^^ dx and k{dx) = k{x) dx := xe"^'^^ dx. (2.27) 

Lemma 2.7 (Green function of killed reflecting Brownian motion with drift). 
The Green function of Z^ is given by 

Gix,y):=\'^^;\^f (2.28) 
\V[x)<p{yj ■■ 0<x<y, 



(2.29) 



where 4),ip : [0,cxi) — ^ M are defined by 

(j){x) = 7re^^Ai(^^ + x), 

i^{x) = Ce'^^'Aiifi^ + x) + e^^Bi(^2 _^ 
with 

C:=CM:=-^^!i4™g. (2.30) 

Proof. As explained in Remark 2.4 the Green function G is obtained in terms of solutions of 

Af^'^^^uix) = 0, a; > 0, (2.31) 

where A'^'^/^ is defined in (2.2G). 

The functions (j) and ip defined in (2.29) are two independent solutions of (2.31). In Lemma 2.8 
we show that they satisfy conditions (i) and (ii) from Remark 2.4, whereas (iii) holds by the choice 
of C = C(/i). It remains to show w{il;,(t)) = 1. Using independence of the Wronskian of x and 
(2.20) we obtain 

ij'iOMO) - ^(O)0'(O) = M^(0)<^(0) + (Bi'ifi^) + GAi'{fi^))nAi{ii^) 

- fi(f>{0)tp{0) - TTAi'{n^){Bi{fi^) + CAiin"^)) 
= vr {Bi'{t?)Ai{n'') - Bi(i?)M'(i?)) 
= Ti{Bi'{Q)Ai{0) - Bi{Q)Ai'{Q)) ^ 1. 
In particular, w{i>, (p) is also independent of /i. Altogether the assertion of the lemma follows. □ 

Lemma 2.8 (Properties of (j) and ?/>). Let and i/j be defined by (2.29). For any > the 

function (j) is strictly decreasing and the function ip is strictly increasing in x. 
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Proof. Properties of the Airy function Ai (see Remark 2.5) imply that for any fi > the function 
(j) is positive and that 0(a;) — )■ as a; — oo. We will show that (p'{x) < 0. To this end, it is 
enough to show that for n > 0, x > 

g{n, x) := ^iAi{^/ + x) + Ai' {^^ +x) <Q. 

First we show that 51 (/i) — g{fJ.,0) < for /i > 0. The assertion is true ior fi = and for 
/i — >■ cx) we have gi{fi) — 0. So ii gi{fj,) is positive on some interval then there is a local maximum 
in some /xq such that on the one hand we have gi{^^o) > and on the other hand 

= g'M = Aii^il) + 2^ilAi'i^il) + 2fioAi"{ij,l) 
= Md) + 2Mo(^^'(Mo) + ^^oAi{pl)) 
^ Ai{fil) + 2filgi{no)>0, 

leading to a contradiction. 

Now we fix /i > and show g{^,x) < for all a: > 0. For a; = it is true by the above 
argument. As x — >■ cxd we have g{fi, x) — > 0. If g{ii, •) has positive values in the interval (0, 00) 
then there is a local maximum a^o such that g{ii, xq) > and 

d 

= 2:0) = fiAi'in^ + xo) + Ai"{n^ + xo) 

= ixAi'{^^ + a;o) + ^i^Ai{^i^ + a;o) + x^Ai{^i^ + xq) 

= ^ig{^JL^ Xo) + x^Ai{^'^ + Xo) > 0, 

leading again to a contradiction. 

It remains to show that ip is increasing. By the choice of C we have ip'{0) — 0. Let h{ij,,x) = 
fiBi{fP + a:) + Bi' {^"^ + a;). For all a; > and /i > we have 

^'(a;) = (/^(m, x) - '^sIm, ^)) > 0. 

5(m,0) 

To see this note that, as we have shown above, g{fi, x)/g{p,,0) — \g{ii, x) / g{fi,0)\ < 1 and 
therefore 

H^J', x) ~ -r^5(M, x) > h{fi, x) - 0) > 0. 

The last inequality follows from the fact that Bi and Bi' are increasing. □ 
We set 

$(a;) := e'^^'Xa;) = 7re-^^Ai(^2 ^ (2.32) 

^'(a;) e-^^^^Al^:) = er"'' {Bi{fi^ + x) + CAi{^i'^ + a;)) (2.33) 

and note that $ and solve the differential equation 

u"{x) + 2iiu'{x) - xu{x) = 0, (2.34) 

and that $ is up to a constant factor the unique decreasing solution of that equation satisfying 
$(a;) — >■ as X — >■ 00. Furthermore we have 

^-'(0) + 2^*(0) '0'(O) = (2.35) 

and simple calculation shows 

0(x)*'(x) - il){x)<^'{x) = w(?A, (j)) = 1. (2.36) 
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Remark 2.9 (Expected killing time and the density of the killing position). 

Using the Green function one can compute the mean killing time of the killed reflecting Brownian 
motion starting in a; > 0. It is given by 

/■oo />oo 

E,[r]=/ G{x,y)m{y)dy= G{x,y)2e-^^'y dy (2.37) 
Jo Jo 

which can be written as 

= 2(c^{x) [ ^{y)dy + ^{x) f <Piy) dy) . (2.38) 



Furthermore the density of Z^_, i.e. the position at killing time is given by (see Borodin and 
Salmincn, 2002, p. 14) 

G{x,y)k{y)^G{x,y)ye''"'y. (2.39) 

□ 

Lemma 2.10 (Exponential moments of the killing position). 

For a < fjL and any x > we have 

Ex[e"^^-] < oo. (2.40) 

Proof. Set y* :— maXj,>o{?/e~'^''~"^^} and recall the function M (for C defined in (2.30)) and its 
bound G* in (2.23) and (2.25). Then (2.40) follows from 



[e"^-- ] = / e"^G(x, y)ye-2^« dy 
Jo 

px poo 

= ct){x) #(y)e-(2A'-")^dy + ^(2;) / y^(y)e-(2M-")yrfy 

Jo Jx 

= e^^7r|Ai(/.t2 + x) y" ye-(^-")^(Si(/i2 + y) + CAi{n^ + y))dy 

+ (Si (^2 ^^^^ CAiifi^ + ?/e-(^-")^Ai(/^2 +y)dy\ 

J X 

< e^^y*M(x) < e'''=y*G*. 



□ 



Lemma 2.11 (Second moments of the killing time). 

There is a positive finite constant C***^'" such that for all x > 

E^[r2] < C'^'ii. (2.41) 

Proof. Since the killing time of the killed reflecting Brownian motion starting in a; > is bounded 
stochastically by the killing time of the killed reflecting Brownian motion starting in 0, we have 

E,[r2] < 1+Eo[t2]. 

By the Kac's moment formula (see e.g. Fitzsimmons and Pitman, 1999, (5) on p. 119) we have 



Eo[t2]=2/ G{0,x)m{x) G{x,y)m{y) dy dx. 
Jo Jo 

Now using again (2.23) and (2.25) as in the proof of Lemma 2.10 we obtain 

/"OO /"OO 

Eo[t^]<AG* G(0,x)m(a;)e^^da; = 8G*7rV'(0) / Aiin"^ + x) dx. 
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Since Ai is decreasing we obtain 



Eo[r2] < 8G*7r7/)(0) / Ai{x) dx 



8G*TT 



where the last equahty follows from (2.21). 

Lemma 2.12 (Bound on the expected killing position starting from x). 
For any a; > 

Proof. By (2.-39) and the definition of $ and ^I* in (2.32) respectively (2.33) we have 



□ 



(2.42) 



E.[Zr-] = / y^G{x,y)e-^^ydy 
Jo 



^^x) 2/2^(2/)e-2^^dy + V(a;) / y^cj){y)e-^^^y dy 

Jo Jx 

px poo 

= (f>{x) y^^{y)dy + i!{x) y^<i>{y)dy. 



Now using the fact that $ and \l/ satisfy (2.34), integration by parts, (2.36) and the fact that 
'i' (f) = jp^ we arrive at 



E,[Zr-] - 0(x) J {yf'iy) + 2/iy«''(y)) dy + i^ix) / (y$"(y) + 2^iy<^'{y)) dy 
= x[4>{x)^' {x) - il:{x)^' {x)] 

-*(a;) + *(0) + 2pLX-^{x) - / *(?;) dy 

ip{x) 



/•oo 

<i>(x) - 2fj,x^{x) -2^ <i>(y) dy 

J X 



a; + 0(a;)*(O) - 2/1 (/)(a;) / '9{y)dy + Tp{x 



$(y) dy 



= X + (l){x)ip{0) - ^j. G{x,y)m{y)dy 
Jo 

= x + (j){x)4;{0)-fiE,[T]. 

Here the next to last equality follows from (2.38) and from V'(O) = *(0)- Since (p is decreasing 
and Ea;[r] < Eq[t] the assertion (2.42) follows. □ 

2.4. Invariant distribution at jump times. In this subsection we consider the increments of 
the Brownian ratchet at the jump times of the boundary process. We show that they constitute 
a Markov chain with unique invariant distribution and compute the expected jump time and the 
expected killing position under the invariant distribution. 

Definition 2.13 (Markov chain at jump times). 

Let {X,TZ) be (7, /z)-Brownian ratchet with sequence of jump times of TZ given by (t„)„>o. We 
define the Markov chain (3^, W, 77) := (F„, Wn, '7n)n=i,2,... of increments at jump times by 

Yn^ Xr^^- Rr„, Wn ^ Rr„ - Rt„_, and ?7"=T«-T„_i. (2.43) 



Since for any k the law of {Y„,Wn,r]„)n=k+i,k+2,... depends on {Yn,Wn,rin)n=i,...,k only through 
Yk, (3^, VV,?/) is indeed a Markov chain. 
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Proposition 2.14. There exists a unique invariant distribution of the Markov chain (3^,yV,77). 

Proof. While uniqueness of an invariant distribution is guaranteed by the coupUng result in 
Lemma 2.3, to prove existence we need to show that the moments of {YmWmrjn) are bounded 
for all n. This implies then tightness of the sequence and also tightness of the Cesaro averages 
of the laws. Weak limits of subsequences of the latter are invariant distributions of the Markov 
chain. Boundedness of the moments of ijn follows from Lemma 2.11 and that of the moments of 
Yn (and Wn since it has the same distribution as 1^) follows inductively from Lemma 2.12. For 
details we refer to the proof of Proposition 5.6 in (Depperschmidt and Pfaffelhuber, 2010). □ 

Proposition 2.15. Let v he the invariant distribution of Yi,Y2, ■ ■ ■ and let denote the ex- 
pectation with respect to that distribution. Then there is a constant K £ (0,oo) so that 

KAYi] = + Ai{ii')). (2.44) 



E.M^"^. (2^45) 

Proof. Let denote the density of v with respect to Lebesgue measure. Let Z^'^ := (z/^'^)t>o 
be killed reflecting Brownian motion with drift starting in random value Y" > with increments 
independent of Y . Furthermore let U be uniformly distributed on (0, 1). Invariance of v implies 
that for killing time t of Z^-^ we have 

As in (Depperschmidt and Pfaffelhuber, 2010, Section 5.3) from that one can obtain the following 
recurrence equation 

/■OO />oo 



f-i^) = /.(2:)y e-'''^G{x,u)du 
and then compute 

fl{z)^-^{z) f,{x)4>{x)dx-<^{z) f,{x)^j{x)dx, (2.46) 

J z Jo 

nOO pZ 

f'l{z) = -^'{z) f,{x)4>{x)dx~<^'{z) f,{x)^{x)dx, (2.47) 

Jz Jo 



f'j'iz) = -^"(z) / fA^mx) dx + v'(z)/.(z)0(z) 



and 



-$(z) / f,Ax)^{x)dx-<^{z)f,{z)i;{z) ^^-^^^ 
Jo 

^~2fif'J{z) + zfUz) + fAz), 
where for the last equality we used (2.34), equations for fl and (2.3G). Thus, 

/:'(z) = -2m/:'(z) + (z/.(z))'. 

Integrating we obtain 

f'J{z) = -2^Ifl{z) + zf^{z). (2.49) 
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The integration constant is zero because from (2.46), (2.47) and (2.35) we see that /"(O) = 
—2fifl{0). By (2.46) the density fi, must be strictly decreasing. Up to a constant factor the 
positive decreasing solution of (2.49) is $ and it follows that 

1 

/^(z) = $(z) with <S>{x)dx. (2.50) 



K 

From (2.49) it follows 

K[Yi]^ xU{x)dx^ I {f:{x) + 2tifl{x))dx 



= -f'M ~ 2Ai/,(0) = -(M»(m') - A^ii?) - 2^,A^{^?)) 

(replace here K by Ktz to get (2.44)) and 

E4iu]=2 Mx) e-^''yG{x,y)dydx 
Jo Jo 

= 2^ f.{x)(4>{x) ^{y)dy + 4>ix) ^{y)dy)dx. 
Now using Fubini's Theorem and then (2.46) we have 

E,[7,i] = 2 / (*(y) / Ux)^{x) dx + $(2;) r U^m^) dx) dy 



2°r/;wl = 2/„(0) = 2M=J'*(^") 







K K 

Again replacing K by Ktt we get to (2.45). □ 



2.5. Regeneration structure. In this subsection we define a regeneration structure of the 
(7, /Lt)-Brownian ratchet and show that the second moments of the regeneration times and of the 
corresponding spatial increments are finite. 

Definition 2.16 (Brownian ratchet as a cumulative process). 

Given a Brownian ratchet [X^IV) = {Xt,Rt)t>o with (Xo,i?o) — {x,0), x > we define a 
sequence of regeneration times as follows 

Po ■■= inf{t >0:Xt^ Rt}, Po ■= inf{i > po ■ Rt- + i?t}, (2.51) 
where Rt- = lmis^t.s<t Rs, and for n > 1 we set 

Pn ■■= mf{i > pn-m : Xt = Rt}, Pn ■■= inf{i > p„ : Rt- ^ Rt}- (2.52) 

Then po < Po < Pi < Pi < • • ■ almost surely and we have po — in the case x — 0. Furthermore 
the sequence (Xp^^-^ — Xp^,pn+i — Pn)n>o is iid. We define 

n 

Mt min{n : p„ > t}, 5„ J^^Xp^ - J and At Xp^ + Xt - Xp^^^ . (2.53) 

1=1 

Then we have 

Xt - Sm, + At, (2.54) 
that is, Xt is a type A cumulative process SMt with remainder At (see Roginsky, 1994). 
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It is well known (see e.g. (Roginsky, 1994; Smith, 1955) cf. also Remark 6.1 in (Depperschmidt 
and PfafFelhuber, 2010)) that to prove the law of large numbers and the central limit theorem 
for {SMt)t>o we need to show that the second moments of pi — po and Xp-^ — Xp^ are bounded 
(this is done in Propositions 2.17 and 2.18). Then, to carry the result over to {Xt)t>o we have to 
prove that the remainder {At)t>o is asymptotically negligible (this is done in Proposition 2.19). 

Proposition 2.17. There exists a positive constant R* such that for all x > 

^ApI] < R* and E,[(pi - pof] < R* . (2.55) 

Proof. In the case re = we have po — 0. Consider the case x > 0. Let Hi be the hitting 
time of of the Brownian motion with drift — /i started in 1. This hitting time has exponential 
moments (we used this fact already in the proof of Lemma 2.3). Furthermore set T ~ 'mi{t > : 
Xt — Rt ^ 1} a-nd note that T < E1/2, where E1/2 is independent exponential random variable 
with rate 1/2, because as long as Xf — Rt > 1 the rate at which the reflection boundary jumps 
into the interval [Xt — l,Xt] (and therefore T occurs) is 1/2. But T also occurs if Xt hits the 
interval [Rt,Rt + !]• It follows that for any initial positions x >0, po is bounded stochastically 
by the sum of independent random variables -E1/2 and Hi, both having exponential moments 
and not depending on x. Thus, E2:[po] is bounded by a constant not depending on x. 

We write Pi — Po = {pi ~ Po) + (Po ~ Po) £^nd argue that each of the terms in the brackets has 
bounded second moments. Note that po — po is the first jump time of the reflection boundary 
after po- Thus, the finiteness of its second moment follows from Lemma 2.11 and that the bound 
there, also does not depend on x. The finiteness of the second moment of [pi — po) follows by 
the same argument as the finiteness of the second moment of po- D 



Proposition 2.18. There exists a positive constant R** < 00 so that for any x >0 

EAXpJ < R** and E^iXp^ - Xp„f] < R** . (2.56) 

Proof. Recall that before touching the reflection boundary the process Xt behaves as a Brownian 
motion with drift — p and is therefore bounded below by and above by a Brownian motion 
without drift. Applying the second Wald identity and Proposition 2.17 (with R* from that 
proposition) we get 

E.[X2J < E,[B%] = E,[po] < ^IMpS < (2.57) 

Now we write Xp-^ — Xp^ = {Xp-^ — Xp^) + {Xp^ — Xp^) and note that as in (2.57) the second 
moment of the first term is bounded by %/ R* . The second moment of the second term is finite 
according to Lemma 2.10 with a bound independent of x. Taking R** to be the larger of these 
two bounds (2.5G) follows. □ 



Proposition 2.19 (Asymptotics of At and Xj^ — i?t). We have 
At „ , X^ - Rt 

—= — > and p > a.s. as t — > 00. 

We omit the proof here since it is almost the same as the proof of Proposition 6.7 in (Dep- 
perschmidt and Pfaffclhubcr, 2010) where the corresponding result was shown for the Brownian 
ratchet without drift. (Note that the definition of ¥„ there should be y„ — supjg[p^_^ p,,]!"''^* ~ 
Xp^ j |. Also the denominator in the last two displays of that proof should be ^/i instead of t.) 
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2.6. Proof of Theorem 1.1. Here we only sketch the proof and refer for details to Section 7 
in (Depperschmidt and Pfaffclhubcr, 2010). 

In the case 7 = 0, the law of large numbers in Theorem 1.1 holds since {Xt)t>o is then a re- 
flecting Brownian motion with negative drift — /i bounded stochastically by a reflecting Brownian 
motion without drift. Therefore Xt/t ^ a.s. as i -> oo. 

Hence, we assume 7 > in the rest of the proof. We use the regeneration structure from 
Definition 2. 10 and set 



r := E.j;[pi - Po], 



X 



(3^ := Var, 



{pi - Po)to 



(2.58) 



Here, r, m and are independent of x due to the regeneration structure. According to Proposi- 
tions 2.17 and 2.18 the temporal and spatial increments pi — po £^nd Xp^ ^ ^po have finite second 
moments. From that and Proposition 2.19 it follows 



Xt t^oo m 

> — a.s. 

t r 



Furthermore, using the CLT for cumulative processes (see e.g. Smith, 1955; Roginsky, 1994) 
and Proposition 2.19 we obtain that for all x e M 

where $ denotes the distribution function of the standard normal distribution, and and r are 
as defined in (2.58). Hence, the central limit theorem holds for a — fijy/r. 

It remains to compute m/r. To this end we use the ratio limit theorem for Harris recurrent 
Markov chains (see e.g. Rcvuz, 1984). Let v denote the invariant distribution for (3^, W,77). 
Using the ratio limit theorem we obtain that 

^=lim^=lim^= lim^= Hm ^=1^ = 

where for the second equality we have used Proposition 2.19. We recall that Ey[Wi] = Ey[Yi]. 

Let V : [0, 00)^ [0, 00), (/i, 7) 1— w(/i, 7), denote the speed m/r as a function of /i and 7. In 
the case 7 = ^ we obtain from Proposition 2.15 

/ 1\ E.[yi] ^ ^lAi{p^)+A^'ip^) ^ If Ai'jp^) X 
"^'"'2; E^[ti] 2A^(/i2) 2\Aiip2)^)' 

Now let /i > and 7 > be given. Using the scaling property (see Lemma 2.2) we obtain 

{^t (^^(2^)2/34 '^(27)2/34 )■ 
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Thus, we have 



^(M:7) = lim 7 



(27)1/3 Hm 



E 



,.1/2,(27)-'/^ 
^(27)2/3* 



(27)i/3.((27)-^/V,^) 

(27)1/3 .^j'((27)-2/3^2 



/■Al' 



2 V Ai((27)-2/3^2) 
which concludes the proof of Theorem 1.1. 



(27)2/3^ 



(27)-^/^) 



3. Ornstein-Uhlenbeck ratchet 

In this section we carry out the same program for the Ornstein-Uhlenbeck ratchet as we did for 
the Brownian ratchet. The arguments in many of the proofs here are similar to the corresponding 
proofs in the previous section. Therefore some proofs in this section will be sketchy. 

3.1. Graphical construction. Let us first recall the definition of the refiecting Ornstein- 
Uhlenbeck (OU) process with infinitesimal drift — /ix and unit variance. Suppose that B — 
(i?(t))oQ is a standard Brownian motion starting in and let xq > 0. Then using the repre- 
sentation of the OU process as a time changed Brownian motion, we obtain that Z := ( Zt 
defined by 

1 



t>o 



Zt := 



-fit 



(3.1) 



is a reflecting OU process with infinitesimal drift —fix and unit variance starting in xq. It is a 
diffusion process on [0,cx)) associated with infinitesimal operator A^^ acting on 



as follows: 



!?(>) :={/eC,2(M+):/'(0+)=0} 



A" fix) ■.^\r{x)~tixf{x). 



(3.2) 



The graphical construction in the following definition (see Figure 3.4) is different from the 
graphical construction of the Brownian ratchet. Here we use a family of independent Brownian 
motions to construct reflected OU processes between the jumps of the ratchet. At any jump 
time a new reflected OU process starts in an initial value chosen uniformly between the state of 
the previous process and zero. Then we stick this "peaces" together to obtain the OU ratchet. 

Definition 3.1 (Graphical construction of the Ornstein-Uhlenbeck ratchet). Let N"' be a Pois- 
son process and let ((i?*^*^ (t))i>o)^_p be independent standard Brownian motions. We de- 
fine a sequence of stopping times (t„)„^o.i,... a sequence of Ornstein-Uhlenbeck processes 
(S^("^)„=o,i,... with = (S't"'')t>T„ reflecting at as follows: 

To = 0, (3.3) 



Si 



(0) _ 



xoe 



—fit 



^5(0)(g2,.t 

2^ 



1) 



Given r„_i and S**^" for some n > 1 we set 

f„ = inf{t > ?„_i : iV' n [0, S^(" 



(3.4) 



(3.5) 
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To 



Tl 



T-3 T4 



Figure 3.4. Graphical construction of the Ornstein-Uhlenbeck ratchet. 



and let z„ be the space component of the almost surely unique element of N'^ f] [0, x {rn} 



Furthermore we set r„ — 



(n-l) 



Zn and for t >Tn we define 



Finally we set for t e [r„,T„+i) 



-Rt = ^ ri and St = 5'^ 



(n) 



and for i > 

Xt = Rt + St- 
Note that by construction wc have 



5 < 



stochastically. 



(3.6) 
(3.7) 

(3.8) 

(3.9) 
□ 



The following lemma is the analogue of Lemma 2.1. Though the graphical construction there 
is somewhat different the proof is similar and will be omitted here. 

Lemma 3.2. The process {X,TZ) := {Xt,Rt)t>Q is a {'j, Ornstein-Uhlenbeck ratchet starting 
in (xq, 0). 

Now we construct a coupling of two Ornstein-Uhlenbeck ratchets starting in ((a;o,0), (a;g,0)), 
where we assume xq > Xq > without loss of generality. Let the Poisson process N"' and a 
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sequence ((Bj )t>o)i=o.i.... of independent standard Brownian motions be given as before. To 
construct the coupling 

({Xt,Rt), {X^,R[))^^^, 

set To = r'o = 0, S'q^'' = x, S'q ■* = x' and define the sequences (tv,)„>o and (Tj'J„>o and the 

^ ^ (n) ~ ^ ^ ^ 

processes (5'''"')„>o and (S' )n>o as in Definition 3.1. Furthermore define S, S' , {Xt,Rt)t>o 
and Rt)t>o as in (3.7). Then define the couphng time by 

Tcoupi ini{t > : St = S'^}. 

Since we can use the same Brownian motions and the same Poisson process for both ratchets we 
have St = S'^ for all t > Tcoupi- Thus, on the event {Tcoupi < oo}, there are n,n' such that for 
fc > 

The following lemma shows that the coupling is successful with probability one. 

Lemma 3.3 (Exponential moments of the coupling time). For any n > there is a > so that 

Proof. We only need to consider the case Xq > x'q. In that case St > S'^ for all t. Furthermore, 
from (3.9) it follows that 5' is stochastically dominated by the reflected Ornstein-Uhlenbeck 
process S^-^K Thus, Tcoupi is stochastically bounded by the hitting time of 0, say Hq, by the 
process S'^^\ For Hq we have (see Borodin and Salniinen, 2002, p. 542) 

¥,^[Ho >t]= Erf(a;o/v/2(e2M* - 1)), 
where Erf (a;) = 2tt-'^/'^ J^^ e""' du = 2n-'^^^ Y.n=o ^^T(^W"- as t ^ oo we obtain 

P.ol^o > t] = ^ +oie-'^') < XoCe-^' 

V"" v2(e2A'* - 1) 

for suitably chosen positive C. 

□ 

3.2. Green function of the killed reflected Ornstein-Uhlenbeck process. In this sub- 
section we carry out analogous computations to those in Subsection 2.3; recall in particular 
Remark 2.4. 

Definition 3.4 (Killed reflecting Ornstein-Uhlenbeck process). Let Z (^t)t>o denote re- 
flecting Ornstein-Uhlenbeck process with infinitesimal drift —fj,x and unit infinitesimal variance 
starting in a; > and let denote the corresponding law on the paths space and the expec- 
tation under this law. Furthermore using an exponentially distributed rate 1 random variable ^ 
independent of Z we define the killing time by 

f = inf{i > : 7 /" Zsds>^}. 

The reflecting Ornstein-Uhlenbeck process killed at rate 'yZ is defined as the process Z^ :— 
(Zf )t>o with Z^ = Zt for t e [0, t) and = A for t > t where A ^ M is the cemetery 
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state. The infinitesimal operator A'^''' corresponding to acts on functions / : [0, oo) M 
satisfying /'(0+) = as follows 

- If'ix) ~ fix fix) - 7x/(x). (3T0) 

The speed measure and the killing measure corresponding to the killed Ornstein-Uhlenbeck 
process are given by 

2 2 

m{dx) — m{x) dx 2e^^^ dx and k{dx) — k{x) dx :— 2jxe^^^ dx. (3.11) 

In the case of Ornstein-Uhlenbeck ratchet the confluent hypergeometric functions play a similar 
role as the Airy functions in the case of Brownian ratchet. In the following remark we collect 
some of their properties that will be needed in the sequel. We refer to (Olver et al., 2010, Ch. 13) 
for most of the properties and for more information on confluent hypergeometric functions. 

Remark 3.5 (Confluent hypergeometric functions). Assume that & ^ Z and consider the Kummer 
equation (also known as the confluent hypergeometric equation) 

xz"{x) + (6 - x)z'{x) ~ az{x) = 0. (3.12) 

Note that solutions also exist in the case b E Z but some of relations that we are going to recall 
in this remark and use in the following may not hold in this case. 

(i) Standard solutions: One standard solution of that equation is given by the function x i— 
M{a, b, x) which is defined by 

n— 

where (a)„ is the Pochhammer's symbol defined by (a)o — 1 and (a)„ = a{a + l) ■ ■ ■ (a + n— 1) for 
n > 1. Another standard solution is given by the function x n- U{a, b, x) which can be defined 
(in the case 6 ^ Z) by 

U{a, b, x) -^Xil^M(a, 6, x) + ^^^~^\ ^-H4{l + a-b,2~b,x), (3.14) 
I (a — + 1) I (a) 

where F denotes the Gamma function. If a 7^ 0, —1, —2, . . . , then U and M are independent 
solutions of (3.12). For a = 0,-1,... we have |F(a)| — 00 and the second summand on the 
right hand side of (3.14) vanishes. Thus, as is easily seen from (3.13) and (3.14), in that case 
both U and M are polynomials which are equal up to a multiplicative constant. The system of 
independent solutions in the case a = 0, — 1, . . . is given by 

U{a,b,z) and z^~^M{a - b + 1,2 - b, z). (3.15) 

(ii) Behaviour in the neighbourhood of zero: For any a, 6 G M 

M{a,b,x) = l + 0{x). (3.16) 

For b G (0, 1) 

U{aAx) = J^^~,^},, +0{x^-'). (3.17) 
i (a — + 1) 

For b G (1,2) 

U{aAx) = + J^'l^^. +Oix'-"). (3.18) 

I (a) I a — + 1) 
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(iii) Behaviour at infinity: For a ^ —1, —2, . . . and & > 

A/(a, 6, x) — — e"^ as a; ^ oo. (3.19) 

r(a) 

Here, as usual, we write g{x) ^ h{x) il g{x)/h{x) — >■ 1 as a; —> oo. Furthermore 

U{a, b, x) ~ x'" as a; cx). (3.20) 

Note also that U is uniquely determined by this property. 

(iv) Bounds for positive zeros of U and M : If a, 6 > then M has no zeros on [0, oo). Let 
P{a,b) be the number of positive zeros of U{a,b,x). First we note that using the Kummer 
transformation 

U{a, b, x) = x^-'''U{a -b+l,2-b,x) (3.21) 

we get 

P{a,b) ^ P{a-b+l,2-b). (3.22) 

If a, b and a + b — 1 are non-integers, 6 < 1 and a + 1 > b then P{a, b) = 0. Let a < and 
1 < 6 < 2. If a;o is a positive zero of U{a,b,x) then by (2.19) in (Gatteschi, 1990) 

xo<M^-a). (3.23) 

(v) Differentiation formulas: We have 

d 



and 



U{a,b,x) ^ -aU{a + l,b+l,x) (3.24) 

dx 

-^M(a,b,x) = yM(a+l,b + l,x). (3.25) 
dx b 

(vi) Recurrence relation: There are many recurrence relations for U and M. We will need 
the following (it follows from (13.4.25) in (Abraniowitz and Stcgun, 1992) and (3.24)) 

U{a, b+l,x)^ U{a, b, x) + aU{a + 1, 5 + 1, a:). (3.26) 

□ 

By straightforward computation one can show that if a function f-y^^{x) is a solution of the 

2 . 

Kummer equation (3.12) with a = —jjjp and b — ^, then 

1 (3.27) 



F,,,ix) = e-^-/'^/,,^ ( + ^x^ 



is a solution of 

= 0. (3.28) 

As in the case of killed Brownian motion with negative drift there are positive solutions (j) and 
ip of (3.28) on [0,oo) with Wronskian w{ip,(p) ip'{x)(f>{x) — il){x)(t>' [x] = 1 such that 

{(/) is decreasing and (t){x) — > as a; cx), ^g-j 
-0 is increasing and '(/''(O) = 0. 
Then the Green function of the killed Ornstein-Uhlenbeck process is given by 

{4,{x)^{y) ■.{)<y<x, 
G{x,y) .= <^ ^ ^ (3.30) 
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We define the function x t-^ p{x) by 

Pix) = -^ + VJ^^- (3-31) 

For general 7 and fi the sohitions of (3T2) with a = — and b — ^ with a = are hard 

to deal with. Thus, we will not formulate an analogue of Lemma 2.7 for the killed reflecting 
Ornstein-Uhlenbeck process. In the following lemma we will identify the decreasing solution (jj. 
For ip we will assume that for any given /i, 7 > there is a suitable linear combination, say M(x), 
of 

M(~i,l,p'ix)] e^nd u(^^,l,p'{x)] (3.32) 

if jg, 7^1,2,..., and 



p{x)m(^^^+^-,^,p'{x)^ and C/(^-^,i,p2(a,)j (3.33) 

2 ^ 
if = 1, 2, . . . , so that the function ip : [0, 00) — > [0, 00) defined by 

^{x) = e~''="/^M(a;) (3.34) 

satisfies the condition in (3.29). 

Lemma 3.6. The function x 1— (f>{x) defined on [0, 00) by 

0(x)=e-^-/''[/(^-^,i,/(x)) (3.35) 
is a positive decreasing solution of (3.2X) with 4i{x) — >■ as x ^ 00. 

Proof. Since (p is of the form (3.27) it is a solution of (3.28). By (3.20) it is clear that 0(x) > 
for large x and (j){x) — > as x — ?> 00. If we show that 4>"{x) > for x > 0, i.e. that </> is convex, 
then it follows that (j) is positive and decreasing on (0, 00). We have 

(j)"{x) = 2^j.x(j)'{x) + 2-fX(j){x) 

= 2^ixe''^2^p{x)U'{~j^, ]^.P^{x)) 



V/\e-^p(x)^C/(l-^,^,/(x)), 



where the last step follows by (3.24). Again by (3.20) it is clear that (l>"{x) > for large x. Thus, 
to show that 4>"{x) > for a; > it is enough to show that this function has no positive zeros. 
Using (3.14) and the definition of M one can easily compute that U{^, |,p^(x)) = and 

U{0, f ,p^(a;)) = 1. Thus in the case ^ € {i, 1} we have 4)"{x) > for a; > 0. 

Let us now consider the case G (0, 1) \ {5}. Then by (3.22) the number of positive zeros 

of X i-> C/(l — |, x) equals the number of positive zeros of a; i-> — j^, ^,x). Since 

1 72 1 
2-4^^ + 1^2' 

by Remark 3.-5(iv), the function x i— > U{1 — |,a;) has no positive zeros, which implies that 
d>"(x) > for x > 0. 
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If > 1 then, by (3.23) all positive zeros of x i->- C/(l — |, x) are bounded by 7^/ — 1. 
Since the set {x > : p^ix) < 7^/^*^ — 1} is empty, the function 0" is positive on (0,oo). □ 

In the following example we compute the function -0 in a special case. 

Example 3.7 {ip{x) in a special case). Assume < < ^, let be as in (3.35) and define 



and 



^(x) ■.= iji{x) + C^ix), with C- ^^^^^ 



0'(O) 

By the choice of C we have tJj'^O) = as required. Furthermore, Remark 3.5(iii) implies that 
ip{x) — >■ 00 for X — >■ 00. To show that 4' satisfies the conditions from (3.29) we need to show 
convexity of -0, i.e. positivity of 'tp" on [0,oo) and i^iO) > 0. 
Using (3.25) we obtain 



^$'{ix) = ^ix$[{x) + 7x^1 (x) = -2/i3/2xe '"p{x)M' (^-^, ^P^i^) 

72 _^ / 7^ 3 2 

-xe I- p{x)M [I- —r,-,p {x) 



(3.36) 



By assumption we have 1 — > and therefore Remark 3.5(iv) implies that the function 
x I— >■ M{1 — 1:?'^(2;)) has no zeros on (0, 00). Since by Remark 3.5(iii) this function is 
positive for large x we obtain ipi{x) > and i^'i{x) > for a; > 0. 

In Lemma 3.6 we have shown that 0" is positive and 0' is negative on (0, 00). Together with 
positivity of 0" and "01 it follows 



0"(x)=<(a;)-ti^0"(x)>O, x>0. 
0'(O) 



Finally we have 



m = ^;,(o) - mm = _mm^^tmm . 

0'(O) 0'(O) 0'(O) 

Here w('0i,0) = '>Pi{x)(j){x) — ^pi{x)(j)' (x) is a constant independent of x, and since tp^lx), 0(x) 
and '>pi{x)4)'{x) are positive for large x and (f>'{x) is negative this constant must be positive. Note 
also that, as is easily computed, w(V'i, 0) = w{'ip, <j)). In order this to be one, we would need to 
normalize by that value. Together with 0'(O) < it follows -0(0) > as required. □ 

Lemma 3.8 (Upper bounds of 4>{x) and ■0(a;)). Let tp and be as defined in (3.34) and (3.35). 
For any positive 7 and fj, there are finite positive Ki and K2 such that 

${x) < ^16'^^/''+'^^' and 0(2;) < K2e-'^'=/''. (3.37) 
2 

Proof. Assume first that 7^ 1, 2, ... . Then two independent solutions of the equation (3.12) 
are given by (3.32). In view of (3.19), (3.20), (3.16) and (3.17) it is clear that any linear 
combination of the functions from (3.32) multiplied by e~'^^^^ is bounded on any interval of the 
form [0, t] for t < 00. Furthermore, as x ^ 00, it grows at most as a constant times e~'^^l ^e^ = 
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g7 /^J?+1x/^^+lJ.x fj.Q]^ -^vhich the first part of the assertion follows. If = 1,2, ... , then two 
independent solutions of the equation (3.12) are given by the functions in (.3.33). Using similar 
arguments as above one can show the first part of the assertion also in this case. 
The second part of the assertion follows from the definition of </> and (3.20). 



We set 

$(x) = e-''^'^(a;), 

and note that $ and solve the differential equation 

1 



-u"{x) = (7a; — n)u{x) — fj,xu'{x). 



From (3.37) it follows 



Furthermore, 



and 



^(x) < Kie"''^' and $(a;) < K2e~'''-'' 
^(a;)$'(a;) - i}{x)%' {x) = (^(x)V^'(x) - ^{x)^'{x) = 1 
$'(0) = 0. 



□ 

(3.38) 
(3.39) 

(3.40) 

(3.41) 
(3.42) 
(3.43) 



Remark 3.9 (Density of the killing position and expected killing time). 

As in Remark 2.9 in the Brownian case, the density of the position at killing time of the killed 
Ornstein-Uhlenbeck process starting in x is given by 



f,{y) = G{x,y)k{y) = G(a:, y)272/e-™ 
The expected killing time is given by 



G{x,y)m{y)dy ^ / G{x,y)2e dy 
Jo 

lix) / ^{y)dy + ii^{x) / ^y)dy 



(3.44) 



(3.45) 



Lemma 3.10 (Exponential moments of the killing position) 
For a < ^/ ji and any x > we have 

E^[e'^^--] < 00. 



(3.46) 



Proof. We have 



/"OO 

E,[e"^--]-/ e"^G(a;,y)27ye-^^«'dy 
Jo 

(x) r ye''y%{y)dy + ii{x) I ye''y^y)dy 



27 



Using (3.37) and (3.41) we see that the term in the brackets is bounded by 



-"fx/p. / ye'yy/t^+°'y dy + e'^^^^^^^^ / ye — 'v/t^^'^y'^+'^y dy 



K1K2 



which itself can be bounded by C(7, ii)xe°'^ for some constant C(7, /i) independent of x. □ 
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Lemma 3.11 (Expected killing position starting from x). 
There is a positive finite constant c such that for all x > 

Ex[Zr-]<x + c. (3.47) 

Proof. Wc have 

E,[Zr-]^2j y'G{x,y)e~^y" dy 

JQ 

pX POO 

= Hx) 2jy^^{y)dy + 4i{x) 2-fy^^y)dy. 



Jo J X 

Let us first consider the two integrals. The functions $ and ^E* satisfy (3.40), which can be 
rewritten as 

27xm(x) = v!\x) + 2fixu'{x) + 2^u{x). (3.48) 
Using this (twice in each computation) and partial integration we obtain 

272/2$(y) dy^ j^y ($"(y) + 2/.y*'(y) + 2y.^{y)^ dy 

= x$'(x) - $(a;) + $(0) + 2fix^^{x) - -$'(x) - —x^{x) 



1 



and 



/ 2-fy^^y)dy^ f y (^"(y) + 2^ly^'{y) + 2^l^y)) dy 

J X J X 



= -x^'{x) + $(x) - 2fix'^${x) + -$'(a;) + — x$(x). 

7 7 

Now from -0$ = 0$, (3.42) and (3.43) it follows 

Ex[Z^-] = {x~-) U{x)^'{x) ~ i}{x)^{x)] + 0(x)$(O) + -0(x)$'(O) 

= x-- + 0(a;)$(O) < X + - + 0(O)$(O) 

7 7 

where the last inequality follows because (f) is decreasing. This concludes the proof. □ 
Lemma 3.12 (Second moment of the killing time). For all x > we have E^It'^] < oo. 
Proof. As in the proof of Lemma 2.f f it is enough to show that Eo[t^] is finite. We have 

/•oo /"OO 

Ea[?^]^2 G{0,x)m{x) G{x,y)m{y) dy dx. (3.49) 



Jo Jo 
It follows that (recall the definition of m in (3. If)) 

1^ r ^ 1 - 

G{x,y)m{y)dy ^ -(j){x) I %lj{y)m{y) dy + -i}{x) \ (l){y)m{y) dy 
^ Jo uj ' 

<-l{x)R\ [ e^y/t" dy+-^^{x)K2 I e"^^/''"^'^' 



rj cy poo 

<-l{x)Ki e^y'^" dy+ -^j{x)K2e"^''f'' e"^^' dy. 

1^ .10 ^ .1 rr 



Now using an estimate for x > 



^lx + + A/tt 2^[i 
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(which can be deduced from (Abramowitz and Stegun, 1992, 7.1T3)), and (3.37) we obtain after 
simple calculations 

G{x,y)m{y)dy < -KiK2{e-^''' ''-{e"'' - 1) + ^e^^/^+'^^'e"''^/''-^^') 

for suitably chosen constant C(7, ji) depending on 7 and /i but independent of x. Putting this 
into equation (3.49) we get 

Eo[r2] <2C2(7,^) 

and the proof is completed. □ 

3.3. Invariant Distribution at jump times. Also for the Ornstein-Uhlenbeck ratchet we 
consider the Markov chain of increments at jump times. 

Definition 3.13 (Markov chain at jump times). Let {X,iZ) be a (7, /i)-Ornstein-Uhlenbeck 
ratchet with sequence of jump times of 7?, given by (t„)„>o. At jump times we define the Markov 
chain (5^,W,77) := (?„, M/,,, ??„)„> 1 by 

Yn = - Rr„, Wn = ^ ^r„„i and ??„ = r„ - f„_i. (3.50) 

Now, using the moment bounds in Lemma 3.3, Lemma 3.11 and Lemma 3.12 the next result 
follows as in the case of the Brownian ratchet. 

Proposition 3.14. There exists a unique invariant distribution of the Markov chain (3^, W,?]). 

Our next aim is to compute the moments under this invariant distribution. First wc show 
that the invariant density of the Y component satisfies a differential equation. 

Lemma 3.15. Let v he the invariant distribution of Yi,Y2, . . . and let be the corresponding 
density. Then f^, is the unique positive decreasing solution of 

i/:'(x) = -fixflix) + jxlix) (3.51) 
satisfying /"(O) = and fv{x) dx — 1. 

Proof. Similar to the case of killed Brownian motion with negative drift one can write down a 
recurrence equation for f^. We have 



and 



fAz)= / fi^ix) / G{x,u)dudx (3.52) 



fliz) = 27(-$(z) j l{x)^{x) dx ~ $(z) fAx^x) dx) (3.53) 

I:{z) = 27(-$'(z) Ux)^{x) dx - $'(z) Ux)i,{x) dx) (3.54) 

C(z) = 27(a;7,(a;))' - 2^x{xfi{x))' . (3.55) 

The equations (3.52)-(3.54) are analogous to the case of the Brownian ratchet with negative 
drift. We only give some details on how we obtain (3.55). 
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Differentiating (3.54) and using the fact that $ and solve (3.40) we obtain 
1 



Jo 

POO 

-2{{jx - fi)^{x) - Mx$'(a;)) / L{z)^{z) dz + ^' {x)l{x)^{x) 

J X 

-2{{^x-^i)^x)-^Jix%'(x)) f U{z)^:{z)dz^^\x)U{x)^{x) 



= t{x){^'{x)^{x)~^'{x)ii{x)) 

+ 2(7x-Ai)(-$(a:) j t{z)^{z) dz - ^x) j t{z)^{z) dz 

+ 2/ix($'(x)^ %{z)^{z)dz + %'{x) U{z)i^{z)dzy 

Thus, using (3.42), (3.53) and (3.54) we obtain 

C(x) = 2jUx) + 2(7x - ^^)fUx) ~ ^tixf^ix) 

= 2j{xflix) + Ux)) - 2fi{xf'Jix) + Jl{x)) 

= 2^{xiix)y -2f,{xfiix)y, 

which shows (3.55). 

Now integrating f'J'{x) we obtain 

£'(x) = 2jxMx) - 2tixfl{x). (3.56) 

Here, the integration constant is zero because by (3.54) and (3.43) it follows that /"(O) = 0. By 
(3.53) X fu{x) is decreasing. This concludes the proof. □ 

Recall p{x) in (3.31). One can check that the general solution of (3.51) is given by t^-t^'-i^ 
multiplied by a linear combination of 

m{-^^^\,\,,^x)) and U[^^^^\\,^x)) (3.57) 
if £;f-i^l,2,...,and 



4^i3 2 



p{x)M[-^ + l,^-,p\x)^ and U {- ^ + ^^-^ {x) ] (3.58) 



if — ^ = 1,2,.... In view of (3.19) the modulus of any solution containing a non-zero 

proportion of M in the above cases behaves (up to a polynomial factor) as e'*'^/'^ as a; — oo and 
therefore cannot converge to for x — > oo. Thus, the density of the invariant distribution of the 
Y component of the Markov chain {y, W, rj) is given by 

= poo , ^ , . , h^^^{x), (3.59) 
Jo ri^,^^{y)dy 
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where 

h,,,{x) := e-^-/^-^-V (^-^ + ^, l,p'(x)^ ■ (3.60) 

The following result shows that defined in (3.59) satisfies conditions stated in Lemma -3. 1-5. 

Lemma 3.16. We have 

(i) hfj_y is positive on [0, oo), 

(tt) /^^^(0)=0, 

(in) /p h^,j{x) dx < oo, 

(iv) h'^ ,^{x) < for X > 0. 

Proof. Let 7,/i, > be given. Throughout the proof we write h for /i^.^. The assertions (ii) and 
(iii) hold because h is a, solution of (3. -51) and because by (3.2(J) h{x) ~ ^-i^/f^-f^^ as x — > 00. 

To show (i) we first observe that by (3.20) h{x) > for sufficiently large x. Thus, it is enough 
to show that h has no positive zeros. By the Kummer transformation (3.21) the function h has 
positive zeros if and only if the function x 1-^ U{1 — , |,p^(x)) has positive zeros which is not 
the case as we have seen in the proof of Lemma 3.6. 

To show (iv) we show that h is convex, i.e. h"{x) > for x > 0. Then (iv) follows from (i) 
and from h{x) — > as a; — )■ 00. For x > we have 



— h"(x) = -tih'(x)+-fh(x) 
2x 

= — /Lt( (— — — 2^x)h{x) 

+ e--/---^2V^(-l, + V7^.)(-(l - ^))Uil - ^, IpHx))) + ^h{x) 
= 2(7 + M^-) (h{x) + e--/---^ (1 - - l^ix))) 

- 2(. + .^.)e-w.-..^ (^(1 _ + (1 _ jL^uil - W 'r''^'^^) 

= 2(7 + M':r)e-^-/^-''-^ - ^,l,p'{x)), 



where the last equality follows from (3.26). Now again by the Kummer transformation h" is 
positive on (0,oo) if and only if the function x 1—^ ?7(— |,p^(x)) is positive on that interval. 
This was shown in Lemma 3.6. □ 



3.4. Proof of Theorem 1.2. The proof in the case of the Ornstein-Uhlenbeck ratchet is almost 
the same as in the case of the Brownian ratchet with negative drift. That is, we can again define 
a sequence of regeneration times and show that the temporal and spatial increments of the 
ratchet between this regeneration times have finite second moments. All the ingredients needed 
for the proof of that have been provided in the previous subsections. We content ourselves with 
computation of the speed of the ratchet. To this end we need (as in Proposition 2.15) to compute 
the expectation of Yi and of fji under the invariant distribution ly. 
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dx 



Using (3. 51) we obtain 

K[Y,] - xl{x) d^=^l°" (f'li^) + 2m:^/:(2^)) 

^^(^-Il{Q)-2^i Ux)dx^=^[-fM-2^,). 

Furthermore, recalhng (3.45), we have using Fubini's Theorem in the second equahty and (3.53) 
in the third 



E.[?yi] =2 Uix) (^0(x) ^iy)dy + ^{x) J Hy)dyj dx 

= 2 / U'(y) / Mx)(j>{x)dx + 'i>{y) / /,(x)i^(x) ) dy 







1 r w ^-(0) 



1 Jq 1 
Now the speed of the Ornstein-Uhlenbeck ratchet is given by 

[fi] _ Ji (0) + 2^Ji _ h'^^^ (0) ir i^) 



^Am] 2/,(o) 2/1^,^(0) /i^^^(o) 

□ 
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